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Введение. В статье изучается осесимметричная задача теории упругости для радиально-неоднородной 
трансверсально-изотропной незамкнутой сферы, не содержащей ни один из полюсов 0 и п. Считается, что 
модули упругости являются линейными функциями от радиуса сферы. Предполагается, что боковая 
поверхность сферы закреплена, а на конических сечениях заданы произвольные напряжения, оставляющие 
сферу в равновесии. 

Целью данной работы является асимптотический анализ задачи теории упругости для радиально-неоднородной 
трансверсально-изотропной сферы малой толщины и исследование на базе этого анализа трехмерного 
напряженно-деформированного состояния. 

Материалы и методы. На основе уравнений теории упругости методом однородных решений и 
асимптотического анализа исследуется трехмерное напряженно-деформированное состояние радиально- 
неоднородной сферы. 

Результаты исследования. После выполнения однородных граничных условий, заданных на боковых 
поверхностях сферы, получено характеристическое уравнение и произведена классификация его корней 
относительно малого параметра, характеризующего толщину сферы. Построены соответствующие 
асимптотические решения, зависящие от корней характеристического уравнения. Показано, что решения, 
соответствующие счетному множеству корней, имеют характер пограничного слоя, локализованного в 
конических срезах. Разветвление корней порождает новые решения, которые характерны только для 
трансверсально-изотропной радиально-неоднородной сферы. Появляется слабозатухающее погранслойное 
решение, которое может проникать глубоко вдали от конических сечений и изменять картину напряженно- 
деформированного состояния. 

Обсуждение и заключения. На основе построенных решений можно определить области применимости 
существующих прикладных теорий и предложить новую более уточненную прикладную теорию для радиально- 
неоднородной трансверсально-изотропной сферической оболочки. 


Ключевые слова: уравнения равновесия, уравнения МЛежандра, радиально-неоднородная сфера, 
характеристические уравнение, погранслойные решения, вариационный принцип, прикладная теория, метод 
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Введение. Одним из свойств материалов, влияющих на напряженно-деформированное состояние 
упругих тел, является их неоднородность. Изучение напряженно-деформированного состояния неоднородных 
тел на основе трехмерных уравнений теории упругости связано со значительными математическими 
трудностями. 

Исследованию трехмерных задач теории упругости для сферы посвящен ряд исследований. 

В работе [1] на основе уравнений теории упругости для сферы получено общее решение, 
удовлетворяющее граничным условиям на контуре в смысле Сен-Венана, проведен анализ напряженно- 
деформированного состояния сферы. В [2] на основе уравнений теории упругости для толстой изотропной 
сферы построены однородные решения, зависящие от корней трансцендентного уравнения. В [3] на основе 
решения трехмерных задач теории упругости для сферы малой толщины изучена точность существующих 
прикладных теорий и дан метод построения уточненных прикладных теорий. В [4] изложена трехмерная 
асимптотическая теория трансверсально-изотропной сферической оболочки малой толщины. В [5] приведен 
анализ трехмерного напряженно-деформированного состояния трехслойной сферы с мягким заполнителем. 
В [6] методом однородных решений изучена задача кручения для радиально-неоднородной трансверсально- 
изотропной сферы малой толщины, когда упругие характеристики меняются линейным, квадратичным и 
обратно квадратичным законами по радиусу. В [7] изучена задача кручения для радиально-слоистой сферы с 
произвольным числом чередующихся жестких и мягких слоев. Показано существование слабозатухающих 
погранслойных решений и возможное нарушение принципа Сен-Венана в его классической формулировке. 
Построена прикладная теория кручения радиально-слоистой сферы, адекватно учитывающая возникающие 
особенности. В [8] с помощью метода конечных элементов и сплайн-коллокации исследована задача теории 
упругости для радиально-неоднородного полого шара. Проведено сравнение результатов, полученных 
методами конечных элементов и сплайн-коллокации. В [9] методом асимптотического интегрирования 
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уравнений теории упругости изучена осесимметричная задача теории упругости для радиально-неоднородной 
трансверсально-изотропной сферы малой толщины. Построены неоднородные и однородные решения. 
Установлен характер напряженно-деформированного состояния. В [10] методом однородных решений 
рассмотрена осесимметричная задача теории упругости для сферы малой толщины с переменными модулями 
упругости. Получены асимптотические формулы для перемещений и напряжений, позволяющие рассчитать 
трехмерное напряженно-деформированное состояние радиально-неоднородной сферы. 

Материалы и методы. Рассматривается деформация в рамках линейной теории упругости 
незамкнутой сферы, материал которой является неоднородным по радиальной координате и трансверсально- 
изотропным. Толщина полой сферы предполагается малой, по сравнению с радиусом и размером по дуговой 
координате. Рассматриваются граничные условия, позволяющие решать задачу в осесимметричной постановке. 
Предполагаем, что сфера не содержит ни один из полюсов 0 и л. В сферической системе координат область, 


занятую сферой, обозначим через Г = {" =["; |], 9е[@; 60,], ФЕ[0; 27}. 
Рассматривается линейная зависимость упругих свойств материала по радиусу: 
(0) ©), — 70) = 
А =ал Г, А, = то Ной ‚ 4 =а> г 2 Аз; = аз Г, Ан =ад Г (1) 


(0) (0) (0) (0) (0) 
где, аи ›а> ›@>> ‚Аз › Ча 


—Щ_ некоторые постоянные величины. 
Система уравнений равновесия при отсутствии массовых сил в сферической системе координат г, 0,9 


имеет вид [11]: 





до, ь, 100% - 20. -0%-бъ +080 





=0, 2 

ОГ г 00 х _ 

00», 100 , 30. +9) 80 (3) 
О’ г 090 т 


где, с 


т? 


Оо Об > Обр —Щ_ компоненты тензора напряжении, которые выражаются через компоненты вектора 


перемещений о, =6, (г,0), о, =0,(т,0) следующим образом [4]: 








о.=А, = + А- ешо +20. + 7 , (5) 
о„=А, — +- Е (бы 5} о. + А›0, се + А,, 7 | (6) 
в = 4 : © +А,)0.+А.0,с1е0+ А, о | (1) 

си у 


Подставляя (5)-(8) в (2)—(3), с учетом (1) получаем уравнения равновесия в перемещениях. 











и ди 
ыр 52 "_ +2 (25% 5%?) =?и, +(250 5% 5% 5% }х 
ди, ди 0°и 0’и, би 
х=? | Очи себ |+ =(Ь®+Ь® ]| - бсюб-+—® +=? РР 40 |= 0, 9 
50 и г) (% и др 8 вер ч | 50? | 20 с (9) 


























д? д ^и ‚ ди 
Бо ы Пе т + (35 + +) 2+ 
др р дро@ 00 
21 (0) - (0) би, (0) (0) 
ч=?Ь® | в срб+ т += ве =? (Боев? 0 + и, =0. (10) 
1 а В и ь 1 7, 5 ы 
Здесь р=-Ш — новый безразмерный переменный; # = т. — | — малый параметр, характеризирующий 
8 п и 
(0) 
з ь 1. о, о _ Ч | 
толщину сферы; х=\ию; ре[-Б Ш и =, и, =—, В, = С безразмерные величины; ©, — 
^ ^ 0 


некоторый параметр, имеющий размерность модуля упругости. 
Предполагаем, что боковая часть границы сферы закреплена, т.е. 
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=0, (11) 
в: (12) 
Считаем, что на торцах сферы (на конических срезах) заданы напряжения 
оо = Л.Р), бы = Л.Р). (13) 
Здесь /,(р), №». (р) (5 =1;2) — достаточно гладкие функции, удовлетворяющие условиям равновесия. 
Решения (9), (10) ищем в виде [3, 4]: 
и,(р,0) =а(р)т(@); и„(р, в) =4(р)т"6), (14) 
где функция т(0) удовлетворяет уравнению Лежандра: 
итеучсшо-и\б)+[ 2 =" @ = 0. (15) 
После подстановки (14) в (9), (10), (11), (12) с учетом (15) получаем: 


Бру зараоунай (вы зы 28)-[2-1Ы? о)- 











(16) 
= - С +50) ар) = (+ +5 250 )а(р)|=0 
д 44 12 44 22 23 12 > 
В атоуеайр- [2 рр ар) р в 
=? (ЗЫ +5 +5 а(р)+=(5% +50) а(р)=0, 
а(р)=0, при р =-1. (18) 
а4(р)=0, при р=+1. (19) 
Решение системы (16), (17) имеет вид: 
ара" [ ре“, + ре С, + рае", + ве С, |, со) 
4(р)=е“ [ле С +4е “*С, чье” С, +4ье “*С я ь (21) 
где С, (п= 1,4) — произвольные постоянные, 
ры? [2-1] + (8 9-48); 
ь = (9 +557), — (255) +52 +53 -№5); 
а, = [5 ЕО = (2% Вр = ); 5; — корни уравнения 
нрырые [2-1 (у +знры вов) вры - ИРЫ «азы? - 
Ра ФО оО 2 + (2 - 1 4 +25 +29 + кр . 62) 


+ в 1 +50 50 +4 -(2 ь ЕН НЫ 0, 


Система линейных алгебраических уравнений относительно С,С,,С.,.С., получается путем 


удовлетворения однородным граничным условиям (18), (19). Равенство нулю определителя этой системы 
является условием существования ненулевых решений и приводит к характеристическому уравнению 
относительно спектрального параметра 2: 


4(2;=) —= (245 — 2.4, (и Р> -БР 57° (=(5 +5) )) + 
+( рр, — р,4, (рн, — рый) 5? (&(5, —5,))=0. 


Уравнение (23) имеет счетное множество корней 2,. Общее решение задачи получается путем 


(23) 


суммирования по корням уравнения (23) 


и, =УМ,а, (рут, (9), 4) 


= 
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=>. м4, (рт) вы 
где р 
брус [реа ето, ре (ета ео} 
а 


—=( +255) =(25.— *), 


+ 4 ы — Р>4! )е 





@ = р>4! (а. ь)е^ 1, (ра Р>4! )е 
©, = Р.& (а, -Ь а“ -Ь (р4> —Р.й )е 


=(я-—25 =(5 +25>) 
› 


г) +4› (Рь — Р.® )е 


—=(25, +55) 


—4, (214> — Р2й де в + Р19› (а, —& )е”* р 
@, =1 (рь — Р›4, мы — Ч (рь = 2в + Рь (а, —& )е = : 


Множество корней уравнения (23) при = —> 0 состоит из счетных множеств корней 


6, 


©, = (24. — Р.4, }е 


=( 251 +5> ) 











д =—и+0(е). (26) 
Для д, имеем: 
1°. При В >0, № -Ь, >0: 
(5—5) (6% 505.5, )з#((я +5,)6) (5, +5) (6% +605,5, эт ((5, -5,)6)=0, (27) 
где 
де; Д-р -Ь, 
в = (25°) в. +00 -5РЫР); В, = ОФ, 
2°. При В >0, Р-Ь, <0: 
выра? — (5 (а? - 8*)-5®) (26а) м 
ва [ 2508? + (65° (а В?) - 5?) | 5п(298)=0, 
где 
я= +62 -Ь, =(а +16); 
ЕЛЬ -Ь, =+(а-18) 
3°. При В >0, №=Ь: 
(5 (256) +2(675° +2) 65=0, (29) 


где 5 = в 
4°. При В <0, №-Ь,>0: 
(5, —5, (5% =, }5й ((5, +5,)6)+ (5 +5, (5% +505, }5й ((5-—5,)56)=0, (30) 
где 
я -Ь; зы Ь,. 
59. При В <0, 6? -Ь, <0: 
выра? (5 (22 - *)- 5 }| зи (26а) + 
(31) 
а 25 в? +(5° (а? °)-5° Е (258)=0, 


где 
х=ы- УР, +а-врх = ее -Ь, а 
60. При В <0, 2 =Ь: 
[В =6%) 51 (256)+2(Ы В) 5? +9) 65 =0 — 2) 
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где $=/Ь| 


Уравнение (27)—(32,) имеет счетное множество решений. 


Приведем асимптотическое построение решений, соответствующих различным группам корней 
характеристического уравнения (23). Подставляя (26) в (24), (25) и, раскладывая полученные выражения по 


степеням & , имеем: 
1 


а) и, (р;0)= УЕ 5 |(ы +5) (5% +5055 5, 605 (6,5, )зт (6,5 р)- 


КЕ 


(5% +552) 5, соз(6,5, )т (6,56) | +0(2) т, (6), 


(33) 


и (р;0)= УЕ‘ 5 (5 +5257 ) (5% +5055 | оз (6,5 )соз(5,5р)-с0$(5,5, )соз (6,552) ] +0(=)} т, (8), (34) 


= 


где д, являются решениями уравнения 


[8 5.) (5 -605 5, (5 +5,)5)+(5. +5, ) (5 +605,5, за (5. —5,)5)=0 


6) и,(р:0) =>. Е (5 +5) (5 +952) зп ( 645, )с0з(дяр)- 


= 


- (5% +5052 5, чп (буя )соз (дур) |+О(=)}т, (0), 


(35) 


(36) 


и, (р;0)= УЕ (52 +05 (5 +5755 )[ (9,5, т (б,5.р) -11 (6,5, 511 (6,5 р)]+0(=}} т, (6), (37) 


К=1 


где бд являются решениями уравнения 


(5, 5.) (6% =), 511 (5, +5)6)-(5 +5, (№ +155 511 (5, —5)6)=0. 
г. 
а) и, (р;0) = УЕ 5 (о +) {[ Вэ (5, 60) 51 (б,ар)- и соз(6,Вр)сй(6,ар)|. 
| авы? соз(5, 6) 51 (6,4) + (5 (в? — В*)-5® (5, В)ей (ба. | = 
-[а зт (5, р) 5й(9,@р)-+ В соз (,Вр)сй(б,ар)] . 


[рарыт чт (5,В)ен(ба) (5% (в? - В°)- 5) сов.) 51 (ба) |+ 0(е)} т(@), 


и, (р;6) = У Е 54 {авы зв (6, вр)св(днар)- (в (а? 1?) 5} 


КЕ 


жсоз(5, 6) 5й (6, @р) | [2авы® с0з(5,В)5й (6,&) +(5> (г -?)-5% 5 (5, В) ей(6,а) | - 


-[2авы? 0$ (68, р)5й (5, ар) ( (“г _ В* ) -5о т (6, Вр)сй (бар) | х 


х[2авы эт (6,в) ей (на) (5 (о? В?) --5Ф со (6,51 (буд) | +0(=)} т, (6), 


где д, являются решениями уравнения 


(38) 


(39) 


(40) 
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в[25тее — (5% (в? - в)- 5%) [58 (26а)+а[ 25? 8" + (5 (=-в?)-5$) (298) =0. (41) 


6) и, (р;0) = УЕ 5 [5 +5) {| В5т(5, Во) сй(6,ор)- и соз( 6, Вр) 51( 6, ар) |. 


| 2аВЫ? соз(6,В)ей(6,а)+ (5? (2 — В?)-5® 5 (6, В)5# (ба) |-- 
—[азт(6, Вр) сй(б,ар) + В соз(5, Вр) эй (б,ар) | 2аВ ВБ” зт(5,В) 51 (6,4) 
-(5 (2 -1°)- 55 )ооз(5, в) св( бла) + 0(ё)|т, (6), 


и, (р:6)=У`Е®5# [[орьо мп(5, Вр) 5й (бар) (5° (а? — в*)- 5х 


КЕ 


хеоз(5,Вр)св(б,ар) | 2аВЫ? соз(5,В)сй(6,а)+(5? (2 — °)--5$ 5 (5,8) 5й (5,2, |- 
-[2авы? со (58, р)сй(9,ар) +(5> (г — В ) -5 )5т(6,6р)5й (бар) | С 


х[2ары? зт(5,В)5й(5,&) =|8 (= _ 6? уе сз (д, В) ей (д.4 | +0(2)} т, (8), 


где д„, являются решениями уравнения 
в] 25 тее — (5? (в? — В°)-55) | #(26а) а 25% в? + (5? (2 -В°)--5$) (298) =0. 
—. 


а) и, (р;0)= УЕ (95? +0) (с03(6,5)с03(6,5р) + рз(6,5) зп(6,5р)) + 


К=1 


(0) ‚2 (0) 
(6 5 5. ) 


50 


т (5,5)с0$ (6,52) + 04) (8), 


2 
(о +59) 


и, (2;0)= > р 
х(р5й1 (5,5) с05(6,5р)- с0$(6,5)з(6,5р))+0(=)} т, (6), 
где бу являются решениями уравнения 
(67 5% (256) +2(6057 +5) 65=0 


В© 52 -5%) 
Е® ( и 44 
6) и, -> и - 


с0$(6,5)т(6,5р) - (55? +50) х 


х(51т (5,5) т (9,52) + рсоз(5,5) соз(5,5р)) + 0(=)} т, (9) 


(7 +) 
(6% +82)65 


и, (р;0) => Е | 


&= 


(42) 


(43) 


(44) 


(45) 


(46) 


(47) 


(48) 
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х| рсоз(6,5)зт (5,5р)-зт(9,5)с05(6,5р) |+ 0(=)} т, (в), (49) 


где д, являются решениями уравнения 
(7-5) т (256)-2(5757 +50) 65 =0 (50) 
4°. В случае В <0, Ь? -Ь, >0 асимитотические формулы для перемещений получаются из (33)-(38) 
заменой 51,5, на 15, ,15.. 
59. В случае Ь < 0, В? -Ь, <0 асимптотические формулы для перемещений получаются (39)(44) 
заменой 51,5, на 15,,15.. 
6°. В случае Ь < 0, 5? =Ь, все асимптотические формулы для перемещений получаются из (45)-(50) 


заменой 5на 15 


Для корней (26) главный член асимптотического решения уравнения (15) при =->0 принимает 
вид [9, 10]: 








- - яч ='\/-52 (0-6, ] (1+0(=)); в окрестности В =8, 
эт 4-92 
т, (в) = ь (© 
ке #288 (9-6, )| (1+0(=)); в окрестности @ =6,. 
тб + 5, . 
Перемещения представим в виде: 
и„(р,0) => Е,а, (рут, (6), (52) 
ри 
и (р, в) = > Е, а, (рт, (в). (53) 
К=1 
Напряжения о и о представим в виде: 
ош = УЕ, (о (рут, (в) +0 (рут (в), (54) 
К=1 
Об в = В о», (р)т, (0), (55) 
К=1 


здесь 
1 ; 1 
ор = [Вау ау [8-1] 


о’ (р) = (5-5) а, (р); 
о 
обр) = 7 [4 Круче ()-4,(Р))] 

Характер решений (33)-(50) существенно зависит от типа корней д. Погранслойные первые члены 
этих решений соответствуют краевому эффекту Сен-Венана [4]. В случае мнимых корней д», эти пограничные 
слои имеют слабое затухание. Таким образом, напряженно-деформированное состояние достаточно далеко от 
торцов существенно от них зависит. То есть в этом случае трансверсально-изотропные свойства неоднородного 
материала значительно, по сравнению изотропным материалом сферы, меняют картину напряженно- 
деформированного состояния. В то же время, при действительных или комплексных д, картина напряженно- 
деформированного состояния неоднородной сферы для таких материалов качественно совпадает, различаясь 
скоростью затухания вышеописанных погранслойных решений Сен-Венана неоднородной плиты. 

Из (51) получается, что при удалении от конических сечений 9=0,(/=1,2) решения (33)-(50) 
экспоненциально убывают. 


Поскольку построенные решения удовлетворяют уравнению равновесия и граничным условиям на 
боковой поверхности, вариационный принцип Лагранжа принимает следующий вид [4, 11]: 


Ле - Л (р))би, +(о„» Л. (Ф))би, | м е’?4р =0. (56) 
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Подставляя (52)-(55) в (56) и считая ОЁ, независимыми вариациями, получаем бесконечную систему 


линейных алгебраических уравнений 
Ува =т» (=12,...), (57) 
К=1 

здесь 


1 2 1 
ан = | 0$ (ра (р)е’*ар| Ут, (6, т! (0, ) |+ | о’ (ра (р)е’арх 
1 5=1 =] 


5 (т, ‘био |+ [оь (ру, пруеар (т (@)т, © 


2 


=>. ти (6,)[ 1. (Ра (р)е?ар-+т,(6,)| №, (рда (р)е??ар | 


5=1 


Система (57) всегда разрешима при физически осмысленных условиях, наложенных на правую 
часть (57). Разрешимость и сходимость метода редукции для (57) доказана в [12]. 

Используя малость параметра = , можно построить асимптотические решения системы (57). 

Результаты исследований. Проведен анализ структуры напряженно-деформированного состояния 
радиально-неоднородной трансверсально- изотропной сферы малой толщины при кинематических условиях на 
боковой поверхности. Показано, что, в случае закрепления боковой поверхности, характер решения 
определяется погранслоями. Оказалось, что асимптотическое разложение напряженного состояния начинается с 
решения описывающего краевой эффект Сен-Венана в теории трансверсально-изотропных неоднородных плит. 
В случае трансверсальной изотропии радиально неоднородного материала сферы некоторые погранслойные 
решения затухают весьма слабо, могут проникать глубоко вдали от конических сечений и изменять картину 
напряженно-деформированного состояния. Выведены асимптотические соотношения для перемещений и 
напряжений, позволяющие рассчитать трехмерное напряженно-деформированное состояние радиально- 
неоднородной трансверсально-изотропной сферы малой толщины с любой наперед заданной точностью. 
Показано, что разветвление корней порождает счетное множество новых решений для трансверсально- 
изотропной радиально-неоднородной сферы. 

Обсуждение и заключения. Асимптотический анализ напряженно-деформированного состояния 
неоднородных оболочек, основанный на трехмерных уравнениях теории упругости, позволяет установить 
границы применения приближенных теорий. Выявленный характер поведения решения вдали от торцов для 
разных граничных условий на боковых поверхностях может стать основой для создания уточненных 
прикладных теорий расчета деформирования радиально-неоднородной  трансверсально-изотропной 
сферической оболочки малой толщины. Одним из приложений проведенного асимптотического анализа может 
служить расчет оболочек с тонкими покрытиями, в которых возникает при этом радиальная 
неоднородность [13, 14]. 
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